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Abstract 

The haaron gas description is reviewed for the QCD vacuum. The role of non- 
renormalizable operators is emphasised in the mechanism which generates 
the string tension. Additional examples are mentioned where certain non- 
renormalizable operators of the bare lagrangian turn out to be important at 
finite energy scale. 



1 Introduction 

The confinement of quarks is a rather mysterious phenomenon of high energy 
physics which challenges our understanding of Quantum Field Theory and 
provides an ever reviving source of inspiration. One may distinguish two 
different confinement mechanisms. The soft one, which is responsible for the 

1 Lecture delivered at the Enrico Fermi International School of Physics on "Selected 
Topics in Non-Perturbative QCD", June 1995, Varenna, Italy. 
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screening of the color charge of an isolated quark by the creation of mesons [|TJ 
in a manner reminiscent of the sparkling of the supercritical vacuum of QED 
around a highly charged ion |2j . This latter occurs when a single electron level 
dives into the Dirac see of negative energy states. The electron of a virtual 
e + e~ pair fills up this hole and the positron escapes to the infinity. The 
semiclassical condition for this to happen is Zcxqed ~ 1- For hadrons Z = 2 
or 3 so one needs a non-perturbative gluonic effect, the anti-screening, to 
raise the running coupling constant to olqcd ~ 1 at the confinement radius. 
This confinement mechanism is called soft because it involves energies at the 
range of the quark rest mass. 

It is widely believed that the gluonic vacuum not only amplifies the cou- 
pling constant at long distances according to the scenario above but produces 
another phenomenon, the hard confinement mechanism. This is the emer- 
gence of the linear string tension which leads to the separation independent 
force between static quark charges. It is a hard mechanism because the 
energy stored in the flux tube is large for well separated charges. 

The two mechanisms compete and it is very difficult to disentangle them 
in the real world which contains the virtual quark-anti quark vacuum polar- 
izations. Only the linear Regge trajectory provides circumstantial evidences 
for the string picture. The more convincing proof of the hard confining mech- 
anism comes from lattice QCD. 

It is clear that both mechanisms are needed to understand the problem 
of quark confinement. But it is the anti-screening and the emergence of the 
string tension in the gluonic vacuum which should be clarified first beacause 
it is simpler. The soft confining mechanism can only be discussed in the 
presence of such a non-perturbative medium. Our attention will be limited 
on the simpler, hard mechanism in the rest of these lectures. 

I think that the well defined and clean environment of the numerical sim- 
ulations of lattice QCD at finite temperature where the deconfinement phase 
transition can be analyzed is the ideal testing ground for our ideas about 
confinement. It is mainly due to the result of such numerical studies that 
the appropriate kinematical framework and the dynamical content of the 
phase transition can be identified. It was found that the center symmetry 
which expresses the invariance of the gluonic system under the fundamental 
group transformations is responsible for the string tension ||. The dynam- 
ical breakdown of this symmetry by the kinetic energy at short time pro- 
cesses links the non-perturbative vacuum with the asymptotically free short 
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distance phenomena and explains some rather unusual features of the high 
temperature phase 0. 

The invariance under the center or the fundamental group transforma- 
tions might be called "top-secret" symmetry after Coleman's beautiful Erice 
lecture f| because it refers not only to the unobservable gauge transforma- 
tions but among them to those which are represented trivially even on the 
gauge field. Only the quark field transforms non-trivially under these trans- 
formations. These transformations are mysterious because it is difficult to 
isolate them in the continuum, can be broken dynamically without loosing the 
global gauge invariance and are not protected by Ward-identities. Another 
characteristic feature of this symmetry is that it requires nontrivial measure 
in the functional integral which is invisible in dimensional regularization. 
Furthermore the typical configurations in the symmetrical realization of the 
path integral contain Dirac delta-type singularities. 

I attempt to reconcile these unusual, sometime confusing features of the 
hard confinement mechanism and perturbative QCD in these lectures. We 
are far from the complete quantitative solution and its outline, the haaron-gas 
vacuum ||, can be given only. The starting point is the gauge invariance, the 
claim that the hard confinement mechanism can be seen only when the gauge 
invariance is guaranteed exactly. In fact, a non-controlled gauge dependent 
component of the vacuum state represents a background charge which can 
screen the well separated test quarks as the virtual quark-anti quark pair 
polarizations do it in the complete vacuum. The ordinary gauge invariance 
is easier to guarantee by satisfying the Ward-identities. The question of the 
center symmetry is more subtle because it is not protected by simple identities 
and can completely be missed in the usual dimensional regularization. Thus 
we shall make sure that the center symmetry is present not only formally 
but dynamically in the low energy effective theory of the vacuum where the 
confining forces can be identified. An interesting question we find in devel- 
oping this effective theory, namely the role of non-renormalizable operators 
will be discussed in some details in the second part of the lectures. 

2 Effective theory of confinement 

This Section is devoted to the isolation and the characterization of the mech- 
anism which is responsible for the linear string tension and the chiral sym- 
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metry breaking in the gluonic vacuum. The center symmetry of the vacuum 
and its consequence is discussed in the first part. A global symmetry does 
not give a detailed enough picture of the dynamics so we need a local effective 
theory for the order parameter. This is introduced in the second part. 

2.1 Center symmetry 

2.1.1 Functional Schrodinger representation 

The center symmetry is simplest to understand in the functional Schrodinger 
representation. The canonical coordinate is the gluon field, A(x) = gA a (x)|?-, 
and the corresponding momentum is the electric field, E(x) = E a (x)^, 

where E a (x) = \ SA t^ ■ The time component of the gluon field is eliminated 
by the choice of the temporal gauge, A (x) = and the hamiltonian is of the 
form 



-Mr/, 



H = -2tr / d 3 x 



r„2 j 



where 



2 V ' 2g 



(1) 



Bidjk = djA k - d k Aj + [Aj , A k ] . (2) 

The gauge fixing and the dynamics are invariant under static gauge trans- 
formations, 

A^x) - Aftx) = iv(x)(d, + A^(x), (3) 

where u>(x) = w(x) G SU(3). We are interested in the gauge invariant 
vacuum sector where the propagator satisfies 

< A f \e- itH \Ai > =< A^|e-^|A, > . (4) 

The usual way to select the gauge invariant contributions from the general 
propagator < Af\e~ ttH \Ai > is to insert a projection operator into the vac- 
uum sector, 

< A f \e- itH \Ai > =< AflPoe-^lAi >= J D[w(x)] < A^|e^|A, >, 

(5) 

where u = e iaaxa and 

V = J D H [a(x)}e i f d3xaa W DEa W. (6) 
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The exponent contains the generator of the gauge transformations, 

DE = <9E + [A, E] (7) 

and the integral variable appears as a static temporal component of the gauge 
field when the path integral representation is worked out for (|5|), « a (x) = 

< A/le-^jAi > = J D H [tgMx)) J D [ A (*, t)]e~ SYM , (8) 

with 

Sym = g^tr J dx(d^A u - d v A„ + A v \f . (9) 

When the projection operator Vq is inserted at each time slice, t n = na, then 
the corresponding integral variable becomes the time dependent temporal 
component of the gauge field, a b (x) = agA^x), 

< A\e- UH \A' > = J D H [agA ( X )} J £>[A(x, t)^. (10) 

We shall see below that the action is given in lattice regularization is this 
case. 

Global gauge transformations act as the basis transformations on the 
gauge field, 

A(x) -> cjA(xV f . (11) 

The center of the gauge group consists of those elements which commute 
with the whole group, it is Zn = {e*~^, i = 1, ■ ■ • , iV} for SU (N). Since the 
center element commute with the generators as well they leave the gauge field 
invariant in (pH|). In SU(2) gauge theory where the center is Z2 = {1,-1} 
and ( |TTD expresses the well known fact that rotations by 2n leave the vectors 
invariant. The relation between the center and the fundamental group can be 
seen by noting that the space of global gauge transformations, SU(N)/Z N , is 
iV-fold connected. The center symmetry is the invariance of the propagator 
< Af\e~ itH I Aj >o under global center transformations, (|j), with cu(x) = e l ~^ . 

It is easy to construct an order parameter for the center symmetry. The 
starting point is to note that the gauge invariance of the path integral ( fL0|) 
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is violated at the initial and the final time slices. In fact, the boundary 
conditions 

A(x,0) =A 4 (x) 

A(x,t) =A/(x) (12) 

do not allow to perform gauge transformations on the time boundaries. The 
remaining symmetries at the initial and the final time slice are the periodic 
gauge transformations with period length t in time since the hamiltonian 
is gauge invariant, [H, DE(x)] = 0. Thus we have no Ward-identities for 
non-periodic gauge transformations. 

The restriction of the space of gauge transformations to the periodic ones 
increases the family of gauge independent variables 0. The observables 
generated in this manner correspond to the Polyakov line, 

n(x) = ?e/» di ' 4(x ' t ' ) , (13) 

which is the path ordered exponential along the stright line connecting iden- 
tical three-space points of the initial and the final time slices. Its eigenval- 
ues are invariant under the gauge transformations which are allowed by the 
boundary conditions, ([12]). The gauge invariant eigenvalues, A a , may serve 
as an order parameter because they transform multiplicatively under center 
transformations, 

A a - e^A a . (14) 

Since the eigenvalues are not distinguishable our order parameter will be 
their sum, 

trft(x) -> e^trft(x). (15) 

It detects the center transformations which are not displayed by the initial 
and final gluon fields, (O). 

The insertion of the projection operator P m t° the amplitude (|j) makes 
the distribution of the order parameter formally center symmetrical. But 
the dynamical preservation of the center symmetry becomes energetically 
less favorable for short time processes. This can be seeby the inspection of 
the effective action for the Polyakov line ||, |5), 

S eff [n] = - In < A^\e- UH \Ai > . (16) 

It is formally center symmetrical but the potential barrier between the center 
symmetrical minima diverges as t — > due to the large kinetic energy needed 



for the finite global rotation of the gauge field in short time [10 . 
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2.1.2 Dimensional v.s. lattice regularization 

The Haar measure is defined by its invariance under group multiplication, 

J d H uf(u) = J d H cuf(u'u). (17) 

This property is required in proving the gauge invariance of @, 

< A w ' | e~ itH | A' > = / D H [u] < A^'- 1 e~ itH \ A' > 
= JD h [uj] < A w |e-^|A' > 

=< A\e- itH \A' > . (18) 

The invariant measure in the path integral (|8|) can be taken into account 
perturbatively. For this end we write the gauge transformation asw = vhv^ 
where h is diagonal, hjk = 6jke mj . The Haar measure reads as 

d H u = dvd N p J 5(£ Uj - 2nn) [J sin 2 Vl^O*. (19) 

n=oo j j<k * 

For the sake of simplicity we shall restrict ourselves to SU(2) gauge theory. 
Then uj = e^^A w jth fi 2 = 1 and 

d H u) = dndpsm 2 ^ = d 3 u— sin 2 ^, (20) 

where dn is the uniform integration over S2. The Haar measure is a rotational 
invariant deformation of the flat integration measure, d 3 u. 

We now return to the path integral where the projection operator is in- 
serted at each time slice and the argument of the Haar measure is agA^x), 

D H [agA b ] = D[n]D[u]e^ 1 d4xlnsin2 au{x)/2 . (21) 
Another form is where the Cartesian coordinate system is kept, 

D H [agA b ] = D[A b ]e^ f d4xln ^? sin2 au{x)/2 

= DK]J A (-^ +0 <" ! )' (22) 

with u 2 = a 2 g 2 A b ) Al is better suited for perturbation expansion. The UV 
divergent factor, is needed in order to make the exponent dimensionless. 
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The nontrivial integral measure can be taken into account by an additional 
potential for the fields in local regularization schemes where the elementary 
volume corresponding to a gauge degree of freedom can be identified. 

The vertices introduced by the integral measure pose a new problem. 
In our effort to implement the center symmetry we introduced a term 0(u 2 ) 
which breaks the gauge invariance of the theory. The remedy for this problem 
actually comes from another 0(u 2 ) piece in the action. As was mentioned 
above the measure term can be treated consistently in lattice regulariza- 
tion only. But the gauge invariance requires the use of the link variables 
Up(x) = e aA ^ instead of the gauge field, A^{x). The expansion of the pla- 
quette in A^(x) gives rise another quadratic piece which cancels the measure 
contributions and restores the gauge invariance at one- loop level. 

Similar cancellations between the lattice vertices continue to occur at 
higher order. A formal proof of renormalizability by induction in the order 



of the loop expansion can be given by the help of the Ward-identities [|11 
There are two kinds of genuine lattice vertices which have no analogy in 
dimensionally regulated perturbative QCD. One, like the measure term, is 
proportional to a negative power of the lattice spacing and UV divergent. 
The other type which is suppressed by a positive power of the lattice spac- 
ing. Neither of these vertices is problematical as far as the overall degree 
of divergence is concerned. In fact, simple power counting shows that there 
are no new types of overall divergences because the dimension of these new 
vertices is supplied by the cut-off itself. To show this consider an observable 
computed in perturbation expansion where we have one coupling constant 
only for simplicity, 

oo 

o = Y,g n in, (23) 

n=0 

where /„ stands for a loop integral. Comparing the mass dimensions of both 
sides we get [O] — n[g\ + [/„]. Since the overall degree of divergence of the 
integral is uj(I) = [I] we have 

u(I n ) = [0]-n[g]. (24) 

This relation shows that the renormalizable coupling constants must have 
nonnegative dimension except when their negative dimension is provided by 
the cut-off itself. The difficult part of the proof is to show that the overlapping 
divergences are removed as well and the counterterms can be chosen to be 
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gauge invariant. The final result is that the contributions of the diverging 
lattice vertices cancel in each finite order of the loop expansion at high energy. 
Furthermore, the order of the loop integration and the removal of the cut-off, 
a — > 0, can be exchanged since the properly regulated theory contains finite, 
uniformly convergent loop integrals. It is worthwhile noting that this holds 
only for theories without anomalies ||12|| , ||14j| . Thus the contributions of the 
lattice vertices cancel and the asymptotically free perturbation expansion 
based on the usual three and four gluon vertices is recovered. 

Dimensional regularization seems to be superior to the lattice regulariza- 
tion because it skips from the very beginning those vertices whose contri- 
butions are ultimately canceled. In fact, the suppressed lattice vertices are 
dropped because only =■ or e-independent pieces show up in the analytical 
regularization. The diverging vertices are absent as well because 

and the right hand side is set to zero. Are the complications of the lattice 
regularization unnecessary ? They are certainly unimportant at high energy 
and in any finite order of the loop expansion. But the situation changes 
when resummation or non-perturbative approximation is sought. The peri- 
odic Haar measure is essential to achieve the vanishing of the Polyakov line 
and confinement at low temperature ||. For any truncation of the Taylor 
expansion for the periodic potential in (|2T|) misses the periodicity and with 
it the confining forces. The perturbative cancellation of the genuine lattice 
vertices is not enough to expel them from the non-perturbative solution. 
But even if we had a non-perturbative argument for the possibility of renor- 
malizing the theory in a gauge invariant manner it would not be enough to 
justify the neglecting of the genuine lattice vertices at finite energy where 
the asymptotic scaling laws do not apply. 

2.1.3 Singular configurations 

Center symmetry implies that the eigenvalues of the Polyakov line, fll3|), are 
distributed equally around the iV-th roots of 1. The local gauge invariance of 
( |IUD allows us to set the integral variables e a9A °^'^ to the unit matrix every- 
where in the space-time except an arbitrarily chosen equal time hypersurface 
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where obviously we have 



e a S Ao(x,t ) = CL ,(x)fi(x)w t (x). (26) 

Whenever an eigenvalue of the Polyakov line is close to a nontrivial iV-th 
root of 1 we must have agAo = O(a ). Thus the fraction of the field 
configurations contain a Dirac delta type singularity, Aq = 0(1/ ag) in the 
continuum limit. In this manner the apparently harmless observation that 
the three vectors remains invariant under rotation by 2tc brings badly sin- 
gular configurations in the renormalized path integral. Since the projection 
operator Vo is inserted at each time slice in (|ToD the configurations dis- 
play the singularities in the fully gauge invariant path integral, too. The 
naive continuum expression for kinetic energy in the Yang-Mills lagrangian, 
— tr(<9oA + DAq) 2 receives the term Y)A from the infinitesimal gauge trans- 
formations performed between the consecutive time slices, 5 A = aDA . This 
expression is not invariant any more when the gauge transformations are in 
a finite, cut-off independent distance from the identity. The more careful 
derivation of the path integral expression which takes such a singular con- 
figurations into account yields the usual lattice regulated kinetic energy in 
terms of the link variable e a5,j4 °( x,i ). The Lorentz invariant extension of the 
lagrangian and the path integral leads unambiguously to Wilson's lattice 
gauge theory. 

The singular structure of the configurations in the path integral is the 
rule rather than an exception. In fact, consider the path integral for a free 
massless particle in D dimension, 

IJ / d0(z)e-3 BD ~ a £W*+^-rt*» 8 . (27) 

The typical configuration is where each contribution to the kinetic energy is 
O(a ), 

0(x + /i) -</>{x) = 0(a^ D/2 ). (28) 

The trajectories are non-differentiable for D = 1, in Quantum Mechanics, 
H|, f[4 |, have finite discontinuities in two dimensions, [|15j, and develop 



Dirac delta type power singularities for D > 2. The Fourier transform of the 
one- loop momentum space UV divergences yields the same result. Due to 
this singular structure the topological concepts introduced on the tree level 



may not survive the renormalization procedure ||T6|. 
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2.2 Effective theory 
2.2.1 Sine-Gordon model 

Our goal is to construct a local effective theory which comprises the low 
energy effects of the center symmetrical fluctuations of the vacuum M. In 
the lattice regulated theory we find the usual three and four gluon vertices, 
infinitely many genuine lattice vertices and the lattice propagator. The lat- 
tice vertices cancel against each other perturbatively in the UV region. We 
suspect that the periodicity of the integral measure is crucial at low energies 
where the cancellation does not hold any more. Thus we shall ignore all but 
the integral measure vertices for A . 

The construction of the effective theory for SU(2) gauge theory is as fol- 
lows: We start in the gauge where A Q is diagonal, Aq(x) = 5 a3 u(x) and a 
U(l) local gauge symmetry is left only. Then we eliminate the off diagonal 
gauge field components which are charged with respect to the diagonal U(l) 
gauge group. The resulting effective theory is a nonlinear U(l) gauge theory. 
So far this is the usual procedure for setting up the Abelian confinement 
scenario, [l7j , except that the Polyakov line is used to single out the Abelian 
subgroup. Then we choose Feynman gauge and eliminate the space-like com- 
ponent of the gauge field. We arrive at an effective theory for u(x) where 
the interactions come from the effective vertices generated by the elimination 
procedure and the measure term. 

Is the measure term is the same in the effective theory as in the original 
bare lagrangian ? Recall that the boundary condition in time for the tra- 
jectories in the path integral fllPf) restricts the gauge invariance for periodic 
gauge transformations. Thus the integral measure for the eigenvalues of the 
Polyakov line is not protected by gauge symmetry and gets renormalized. 
But in the same time the center symmetry remains present and requires that 
the integral measure be periodic. So the integral measure for the diagonal 
component of the local Ao(x, t) must get renormalized as well. In fact, the 
eigenvalues of the Polyakov line are given by the field u(x), 

A (x) = e ± ^/* u ^, (29) 

and its integral measure would not be renormalized unless the integral mea- 
sure for u(x) is renormalized. 

In the spirit of the gradient expansion we ignore all derivative coupling 
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and truncate the effective theory onto the lagrangian 

L= l -{d„u) 2 -V{u) (30) 
where the local potential is periodic, 

V(u + 2tt/k) =V(u), (31) 
V(u) = 2 ^2 v m cosrriK,u. (32) 

m>0 

The contribution m = is eliminated by the condition / duV(u) = 0. This 
four dimensional sine-Gordon type model is our effective theory. It will be 
shown that despite the absence of the other interaction it reproduces the 
salient non-perturbative features of the vacuum. 

2.2.2 Haaron gas 

Consider first a simplified version of the effective theory with one Fourier 
mode only, 

L = -{d^uf — 2Acos/t-u. (33) 
The generator functional defined as 

what is expanded in A, 

oo \n n . 
n=0 j=l J m(j)=±l 

y , ^[ w ] e _ / d4x [i^ u W) 2+i9j W u Wl +iK SJ=i m O') u ( :i; j) 

oo yn n . 

= J2 — f II ( / <^ X 3 S ) I mmIjjImJ, (35) 
n=0 U - 1=1 J m(i)=±l 



where 



T - p _ tEI , =1 m 0')G'(a; J -x fc )m(fc) 

r -p-4-J d 4 xd 4 yJ(x)G(x-y)J(y) 
1 J J ° j 

j mJ _ e -9 K Z)"=i m 0')/ d4 f G ( :E J~f) J (f) (3g) 
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and 

d 4 p ptpx 



G{X) = J (2^^ = (3?) 

The interpretation of these expressions is quite straightforward: The sine- 
Gordon model is equivalent with the grand canonical ensemble of point par- 
ticles with charge ±k and fugacity A interacting via the Coulomb potential 
fl3"7|). The logarithm of the factors I mm , Ijj and I m j contains the energies due 
to the self interaction of the Coulomb particles, the source and the energy of 
the source in the Coulomb fields of the particles, respectively. 

We now return to the effective theory ( p0[) where the repetition the steps 
shown above gives 

Z[J\ = J I)[M]e _ ^ d4l '^ eMU ' ;t '' 2+ ' 9j ' 3: ' u ' 3; '' +2 ^™> ''"' cosmK ' 1 

oo n „ oo 

= S~II( \ <& x 3 J2 v m(j))ImmIjjImJ- (38) 
n=0 n - j=l J m(j)=-oo 

The charge is quantized in the units of k and the Fourier coefficients turn 
into the fugacities. The Coulomb particles are called haarons because they 
represent the effects of the Haar measure. 

The haaron-haaron interactions polarize the gas and introduce the Thomas- 
Fermi screening. The screening mass square is the curvature of the effective 
potential of the field u{x) at the vacuum expectation value. The mani- 
fest center symmetry, the invariance of the theory with respect to the shift 
u(x) — > u(x) +2tt/k, allows the constant effective potential only. This is in 
agreement with earlier numerical finding, namely that the eigenvalues of the 
Polyakov line decouple from the dynamics in the low temperature phase of 
QCD and its distribution is given solely by the integration measure without 



the influence of the action [18[]. The absence of the screening in a Coulomb 
gas is attributed to the negative fugacities. Some of the Fourier coefficients, 
v m , can be negative indicating negative probability in the classical Coulomb 
gas picture of the sine-Gordon model. This makes the cancellation between 
the screening effects of different charges possible. In other words, the self 
energy computed in the framework of the perturbation expansion may be 
vanishing for appropriately chosen coupling constants, v m . The dynamical 
issue of confinement is the assumption of the center symmetry, the absence 
of screening in the effective theory. By the help of the massless propagator, 
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fl37|), the confining forces and the chiral symmetry breaking can be obtained 
in a partial resummation of the perturbation expansion. 

After having excluded the generation of the screening mass the polar- 
ization of the Coulomb gas becomes less important and the approximation 
Imm = 1, the neglect of the Coulomb interactions between the particles is 
more reasonable. The remaining contributions can be resummed in (|38|) 
yielding 

Z[J] =e~^$ d 4 xd i yJ(x)G(x-y)J(y)- J d*xV(ig J d 4 yG(x-y) J(y)) 

= e ~ J d 4 x(p(x)U(x)+V(iU(x)) _ ( 3 g) 

The exponent of the second equation gives the energy of the source system 
as the sum of the perturbative current- current interaction energy and the 
potential energy in terms of the Coulomb potential created by the source, 

U(x)=gJ d*rG(x-y)J(y). (40) 

The only unusual detail is the factor i in the argument of the potential V 
whose role is to remove the periodicity after the elimination of the quantum 
fluctuations. It is worthwhile noting that the Wick rotation into real time is 
nontrivial and the factor i remains present. 



2.2.3 Static charges 

The energy E of a test charge can be read off in the long time limit as 
tE = — In Z by choosing 

J{x)= Jds^ 4 \x-y(s)), (41) 

where ?/ M (s) is the world line of the charge. We take a static charge, 2/ M (s) = 
(s, x ), and use ([39]). 

for the non-perturbative second term in the exponent of (|39|). 

This result is typical inasmuch as the argument of the potential V(u) is 
small in the IR since the Coulomb field created by the source approaches 
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zero at the infinity for localized source. Thus the infrared contribution to 
the energy is controlled by the behavior of the potential around zero, V(u) = 
V(0) + ^V"(0) + 0{u A ). Hence we have 

E np = -V"(0) (£) 2 J d 3 x± + const. (43) 

The curvature of the bare potential should be negative, V"(0) < 0, in order 
to arrive at flat effective potential in the IR. The linear infrared divergence 
of the integral indicates the absence of localized charged among the asymp- 
totic states of the effective theory. Note that the Thomas-Fermi screening 
would regulate this integral and localized charges would be allowed. Thus 
one expects the effective theory to undergo a phase transition and develop 
screening at the deconfining phase transition. 

Similar computation allows us to derive the static potential between test 
charges. To this end we take 



J(x) = J ds[6^(x - y{s)) - 6 (4 \x - y(s) - L)], (44) 
with L = (0,L). The non-perturbative part of ( p9|) gives 



E np = - I d 6 xV 



ign ign 



-nO)yL| + $>|L|-n (45) 



v 47t|x| 47t|x — L| 

<r 

-|-Lf| -f 

n=0 

The leading infrared part is a linearly rising potential with the string tension 

a = -V"(0)^, (46) 
and the coefficients in the sub- leading pieces, c n , are ultraviolet divergent. 

2.2.4 Nambu-Jona-Lasinio model 

It is straightforward to include quark fields into the effective theory. The 
resulting lagrangian is 

L = l(d^u) 2 - V(u) - #<? - S7ocA#. (47) 
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It is not Lorentz or gauge invariant since the diagonal component of Aq are 
kept only after making approximations. The elimination of the field u{x) in 
the free haaron gas approximation gives except the external source is 
replaced by J(x) = ip r yo<T 3 ip, 



S = J d 4 x j- W + V(ig J d*yG(x - y)J{y)) 

2 

- 9 - J d 4 xd 4 yJ(x)G(x-y)J(y) 

d 4 x{-#^ + V{iU{x)) - | j d 4 xJ(x)U(x)y (48) 

where (fiTf ) is used to express the Coulomb field of the dynamical source J. 
The second term in the first equation gives an infinite series of non-local 
vertices. The 0(J 2 ) action is 

S = - J d 4 x^tfj 

+ J d^yJ{x)\- 9 ^G{x - y) + ^G 2 (x - y)}j(y), (49) 

where 

G ^/(W- < 50 > 

Notice that the static potential generated by G 2 (x) is linear with string 
tension a. There have been several studies of the NJL model with such kind 
of interaction and the self consistent gap equation approximation shows the 
dynamical breakdown of the chiral symmetry induced by the strong repulsion 
[0- 

It is illuminating to eliminate the quark field degrees of freedom in favor 
of a composite meson field. By the help of the identity 

D[&]e~H dxd V^)K' 1 (^,y)Hy)+J dxA(x)®{x) __ & J dxdyA(x)K(x,y)A(y) 

we find 

S = -tr m[i<? + 7o a 3 $] + — J dx^^^^. (52) 
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Notice the strong IR dependence of the wavefunction renormalization con- 
stant 

Z(P 2 ) = ^J— 2 - (53) 

Since the restoring force for the small fluctuations of the meson field is weak 
in the IR there are strong interactions between the low energy mesons and 
the quark-anti quark vacuum polarizations. 

It is worthwhile noting that the Euclidean mesonic effective theory is 
defined only below the Landau pole p 2 < ^fcr since beyond this limit the 
Hubbard-Stratanovich transformation (|5T| ) requires imaginary meson field. 

2.2.5 Haaron gas and localization 

There is a simple physical picture behind the haaron gas description of the 
QCD vacuum. Consider the quenched multi-quark Green function in the free 
haaron gas approximation, 

< 0| T${ yi ) ■ ■ ■^{y n )^{z 1 ) ■ • -#Zn)]|0 > 

oo j n p oo 

= ErjIl(MS E v mU)) 

n=0 U - j=l J m(j)=-oo 

< ^{yr) ■ ■ -4){y n )ij){zx) ■ ■ "4}(z n ) > U} (54) 

where the factor in the last line is the free fermion Green function on the 
imaginary background field 

fM = M7^I. (55) 

e.g. 

< $(y)iP(z) >u= \- 1 —— -] (y, z). (56) 

The expression (|54]) is reminiscent of the fermionic Green functions in 
the dilute instanton gas approximation. But our result is not a semiclassical 
contribution, the haaron gas may be dense and is not an extremum of the 
Yang-Mills action. Furthermore there are no zero modes and result is infrared 
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finite. For the dilute haaron gas the fermion propagator can be factorized 
and the zero mode dominance yields the condensate 



< 0|^|0 >= -V(0)tr / d A x < V>(aOV>(z) >u (57) 



where propagator in the right hand side is evaluated on a single haaron 
background. Similarly to the case of the instantons there are localized zero 
modes for the Dirac operator which make the integral finite and generate the 
chiral symmetry breaking. 

The hard confinement mechanism is similar to localization observed in 
strongly correlated electron systems, [pOf| , except that it takes place in the 



space-time rather than the three-space. We have seen that the IR divergence 
of the single quark energy makes the quark propagator vanishing. It is easy 
to reproduce this result by the localization scenario. In fact, consider an 
isolated quark which is sent through the vacuum. The oscillation due to the 
phase shift which is generated by the long range haaron potential cancels the 
propagator. 

In order to understand the propagation of a quark-anti quark pair we 
have to take into account the correlation between the haarons. The coupling 
constants v m provide the scale for the finite haaron density in the vacuum and 
the average haaron distance appears as a correlation length of the quenched 
Coulomb potential generated by the haarons. Let us follow a very simple 
way of keeping track of this correlation: The haaron potential is considered 
completely correlated or uncorrelated for separations which are less or more 
than the correlation length, respectively. The propagator of a meson can be 
written as the sum over the world lines of the quark-anti quark pairs. The 
phase shifts of the quark-anti quark state cancel so long as the separation of 
the quark-anti quark pair is smaller than the correlation length. When the 
distance between the pair increases beyond the correlation length then the 
statistically uncorrelated potential creates a non-vanishing phase shift which 
in turn suppresses the contribution after the haaron gas averaging. Thus the 
quark-anti quark pair tends to stick together and the confinement radius is 
the correlation length of the haaron potential. The similarity between this 
scenario and the realization of confinement in the stochastic vacuum [^TJ is 
remarkable. 
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2.2.6 Double role of V(u) 

The periodic potential V(u) has two rather different role in our approximate 
solution. On the one hand, it controls the large amplitude fluctuations of the 
field u(x) in the effective lagrangian, (^). The order of magnitude of the 
center symmetrical fluctuations is - and the behavior of V(u) in this range, 
in particular the periodicity, is important. On the other hand, the same 
potential appears after the partial resummation in the free haaron gas. Its 
the typical argument is purely imaginary and small in absolute magnitude as 
far as the IR physics is concerned because the Coulomb potential approaches 
zero for large distances. The approximation V(u) ~ V(0)+V"(0)u 2 /2 around 
the origin is sufficient to express the string tension and the chiral condensate. 
After the elimination of the fluctuations the non-fluctuating argument of the 
potential is small and the IR physics is governed by the behavior of the 
potential around zero. 

One finds that the periodicity of the potential and the manifest center 
symmetry has two consequences in the effective theory: In the UV it keeps the 
fluctuations of the field u(x) large and thereby it creates the non-perturbative 
environment. In the IR regime it protects against mass generation. This lat- 
ter appears as the discrete analog of the chiral symmetry. What is further- 
more interesting is that a discrete symmetry is responsible for the massless 
behavior. 

2.2.7 Center symmetry and instantons 

The perturbative ordered vacuum is based on the configuration = and 
the small fluctuations around it. This gives y / <i 3 xtrfi(x) « 2 in SU(2) the- 
ory which indicates that the transition amplitudes are not center symmetri- 
cal. What modes are responsible for the restoration of the center symmetrical 
vacuum, y J d 3 xtifl(x.) = 0? It is reasonable to expect that localized config- 
urations which have large entropy and interpolate between the center related 
minima of the effective potential for the order parameter will reach this goal 
for sufficiently large values of the time parameter t of the matrix element of 
the time evolution operator (|10|). 

The shape of such an interpolating configuration can be found as the 
solution of the Yang-Mills equations of motion. This turns out to be an 
instanton. The reason is that the map of the three-space into the gauge 
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group is S3 — > S3 when the space is spherically compactified. The Polyakov 
line configurations are labeled by the Pontryagin index, the winding num- 



ber. But this is just the topological charge, ]2"2"| . Thus the Polyakov line of 
an instanton winds around the SU(2) group space as the space coordinates 
moves around the whole three-space. The asymptotical value for large co- 
ordinates is Q ~ 1 and in order to cover the whole group space it has to 
take the value Q — — 1 somewhere. Due to the rotational symmetry of the 
solution this happens at the center of the instanton. In this manner trfi(x) is 
a spherically symmetrical localized function which interpolates between ±2. 
The gas of such "domains" restores the center symmetry. Note that the size 
of the instantons which might play role in the center symmetry restoration 
must be around the confinement radius. 



3 Universality and condensates 

The picture of the confining gluonic vacuum outlined above raises more ques- 
tions than answers. The justification of the emphasis put on the center sym- 
metry comes from the numerical experiences in lattice gauge theory, namely 
from the observation that the fate of the center symmetry is related to the 
existence of the string tension. But the important IR effects are produced 
by the Haar measure vertices of the path integral, by the vertices which are 
labeled as non-renormalizable or irrelevant according to the renormalization 
group. These vertices are set to zero in dimensional regularization. Does 
that mean that the usual class of renormalizable field theories is not suffi- 
cient to parametrize all possible physics of a given set of particles and there 
might be different QCDs? We are not in the position to answer affirmatively 
this question which goes beyond the perturbation expansion. Instead, my 
goal will be to indicate only a gap in the usual argument which might make 
possible to describe the non-perturbative IR effects more systematically and 
in the same time to find new continuum field theories. 

There is a well known example where the usual power counting might 
miss a relevant operator, the strong coupling QED p3 |. The point is that 



the anomalous dimensions should be included into the usual power counting 
argument in classifying the operators of the theory. The extrapolation of the 
perturbative anomalous dimension of the electron field gives dimension four 
for the four fermion operator when e ~ 1. Thus we have the possibility for a 
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new relevant coupling constant in strong coupling QED. 

We shall give two other examples for the generation of new relevant cou- 
pling constants. They might be more realistic because one of them refers to 
asymptotically free theories and the other is the Higgs sector of the Stan- 
dard Model. In both cases the apparent violation of the universality and the 
emergence of new parameters are related to condensates. 



3.1 Localized saddle points 

Our example for non-perturbatively generated relevant coupling constants in 
asymptotically free models is based on the higher order derivative terms in the 
lagrangian. If localized saddle points, coherent states, appears in the theory 
then the higher order derivatives may deform them substantially. Depending 
on the sign of the derivatives the saddle points may shrink to the cut-off size 
and saturate the path integral with cut-off effects at each length scale. What 
is interesting is that this may happen despite the weakness of the coupling 
constant at the cut-off scale because the derivatives of the saddle point are 
sufficiently large. 

3.1.1 Bare vs. renormalized expansion 

The starting point is the difference between the bare and the renormalized 
saddle point expansion [p4j| . The well defined path integral is given only for 



the bare, regulated theory, 

= J D[A»}e- s °\l -S l -S C T + -- j, (58) 

where the small parameter of the perturbation expansion is The bare ac- 
tion is the sum of the renormalized one, Sr = So + Si, and the counterterms, 
Sct- The new expansion parameter, g^, is obtained in the renormalized per- 
turbation expansion after taking into account the cancellations in the Taylor 
expansion of the second equation. But note that this procedure is formal in 
the sense that there is no well defined path integral with the small expansion 
parameter g^. In such a manner the reliability and the applicability of the 
perturbation expansion should be investigated in the framework of the bare 
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rather than the renormalized series since the expansion is done before the 
cancellations take place. In fact, the phase transitions provide evidences that 
the renormalizable running coupling constants can not even characterize the 
theory in a unique manner. 

The difference between these two expansion schemes is more pronounced 
for theories with dimensionless coupling constants only, such as the SU(2) 
Yang-Mills theory. Here the renormalized saddle point expansion [25|| is based 
on the strategy that the saddle points are selected by Sr and the countert- 
erms are taken into account on the higher loop level only. The tree level 
saddle points are the instantons whose action is independent of the instan- 
ton size due to the scale invariance of Sr. The higher loop contributions 
break this scale invariance in a manner that large instantons are preferred 
and the non-interacting instanton gas picture becomes inconsistent. The bare 
saddle point expansion brakes the scale invariance on the tree level already 
since a scale parameter, the cut-off, appears explicitly in the regulated action. 
Such a tree level violation of the scale invariance is governed by the dimen- 
sional parameters of the bare action. In order to explore the possibilities of 
the breakdown of the scale invariance we introduce additional dimensional 
coupling constants in the bare theory. The coupling constants with positive 
mass dimension are superrenormalizable and influence the distribution of the 
large, i.e. cut-off independent saddle points. They are excluded in Yang-Mills 
theory by symmetry. The coupling constants with negative mass dimension 
are non-renormalizable and govern the distribution of the ultraviolet saddle 
points around the scale of the cut-off. 

The regulators are represented by irrelevant operators because they are 
supposed to supress the fluctuaitons in the UV regime only. In this manner 
different regulators give different distribution for the small instantons. The 
instantons whose size parameter is in the vicinity of the cut-off and whose 
dynamics is influenced by the regulators will be called mini-instantons. The 
usual Wilson-type single plaquette action decreases monotonically with the 
size of the instantons. The too low action of the mini-instantons leads to di- 
vergent topological susceptibility and nonunique topological charge in lattice 
regularization. The usual strategy to cope with this problem is to construct 
topological charge operator or improved action which cuts out or suppresses 
the mini-instantons, alas topological defects This is certainly justified 
when we are to establish a lattice regulated Yang-Mills model as close as 
possible to the continuum perturbation expansion. But our goal is different 
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in the search of the confinement mechanism when the results of the dimen- 
sional regulated renormalized saddle point expansion can not always be used. 
Instead, we take the bare theory seriously, on a non-perturbative manner at 
each length scale and try to trace downs the modifications of the dynamics 
due to the mini-instantons, if exist, even in the IR regime. 

3.1.2 Mini-instantons 

We start with the bare theory which involves some higher dimensional terms, 
L B = -±- 2 F; u (5 ab + ^D 2ab + 2L D ^ F w. (59) 

The higher order pieces usually come from a heavy particle exchange in the 
energy range A. The coupling constants Cq/ £1X6 irrelevant according to the 
power counting. The simplest way to see this is to consider an observable 
obtained in the perturbation expansion, 

<o>=^ 0( ^c4,c4f)=^E^(^) l (^)'wf). 

jk£ 

(60) 

The characteristic scale, /x, of the observable is used to give its dimension and 
each insertion of a vertex with the coupling constant jfe brings the factor fi a 
by dimensional reasons. So long as the theory is infrared finite and the limit 
^ — > is convergent the c a dependence drops in the renormalized observables. 
The power counting gives another important result, which has already been 
mentioned above in connection with the renormalizability of lattice gauge 
theory. Namely, the coupling constants c a do not harm renormalizability, 
they are actually a variant of the Pauli-Villars regulators. 

The power counting argument is no longer valid if there are localized 
saddle points in the theory. In that case the dependence in the coupling 
constants is not necesseraly polynomial and the suppressing factor (^) a may 
be missing. As an example take Zi/Zq, the ratio of the partition function 
of the one and the zero instanton sector. In order to find the saddle point 
consider a one parameter family of the instanton configurations labeled by a 
scale parameter p. The action of the instanton is 

s ^ = -^-T^ + i^) (61) 
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where c a /c a , a = 2,4 are positive, p and g independent constants. The 
stable mini-instanton size, p, is obtained by solving (/5A) 2 = 204/62, 

SbW = %(l-£\ (62) 



g 2 V 45 4 

The one loop approximation yields 

Zi 



„ Cg p VA*e- SB ®, (63) 
^0 



where V is the four volume ||27|| . Observe the absence of the perturbative 
suppressing factor l/(Ap) 2 in ( [32]) and (|63[) . What is not suppressed here is 
a cut-off contribution because l/(Ap) 2 = 0(A°). 

A more careful analysis of the fluctuation determinant shows that the 
mini-instantons dominate (|63|) for 1 — < X < 1, where X is a given 
positive constant. When this inequality is satisfied then one of the following 
two possibilities is realized: (i) c a are relevant or (ii) c a are irrelevant but 
the beta function for g is non-universal. In fact, let us assume that c a are 



irrelevant. Then the cut-off independence of (33) gives the beta function 



In either case the coupling constants c a modify the physical content of the 
theory at length scales which are independent of the cut-off. In other words, 
the saturation of the path integral by the mini-instantons which are close 
to the cut-off makes the cut-off effects present and changes the dynamics at 
finite length scales. 



3.2 Multiple fixed points 

Another example for the unusual relevant parameters is when there are sev- 
eral fixed points in a theory and a given operator is relevant at one fixed 
point and irrelevant at another one. The appearance of multiple fixed points 
is typical in Particle Physics where different interactions are found to be 
dominant at different energy scales. 
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Statistical Physics Quantum Field Theory 

U.V. fixed point Renormalized theory 

Irrelevant coupling constant Non-renormalizable coupling constant 

Relevant coupling constant Renormalizable coupling constant 

Universality Considering renormalizable theories only 



3.2.1 Theory of Everything 

Let us imagine the renormalized trajectory of the Theory of Everything 
(TOE). This theory includes all physics, at any length scale, by definition. In 
order to understand its features it is useful to summarize the connection be- 
tween the language of the renormalization group when applied in Statistical 
Physics and in Quantum Field Theory. 

The UV fixed point is where the correlation length is infinite. It cor- 
responds to the infinite value of the cut-off, to the renormalized theories. 
The applicability of the linearized version of the blocking relations gives the 
scaling regimes. In these regions we have a classification of the coupling 
constants. Those coupling constants which decrease or increase when we 
move towards the IR regime are called irrelevant or relevant, respectively. 
The irrelevant coupling constants increase as we move into the opposite, UV, 
direction. This is what we always do in the renormalization of a Quantum 
Field Theory. Thus the irrelevant operators prevent us from going "back" 
to the UV fixed point and from removing the cut-off. Hence they are non- 
renormalizable. 

Finally, universality states that the physics at the length scale which 
corresponds to the IR end of the scaling regime is insensitive for the choice of 
the irrelevant coupling constants. This translates into the usual rule of model 
building in Particle Physics, the use of renormalizable theories only. The non- 
renormalizable theories were first excluded due to their uncontrollable UV 
behavior. By some tricks, such as the use of the cut-off to suppress the non- 
renromalizable coupling constants mentioned above may help to eliminate 
the divergences from the theory but the predictive power and the simplicity 
are lost. A simpler argument to ignore non-renormalizable coupling constants 
is offered by the renormalization group: They become small anyhow at the 
scale of the observations so we might as well set all of them zero from the 
very beginning, at the cut-off. 

We now return to the renormalized trajectory of the TOE. At its asymp- 
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totically high energy scaling regime, in the vicinity of the UV fixed point the 
trajectory is governed by the relevant coupling constants. But it is the lesson 
of the Wilson-Kadanoff blocking that all coupling constants what might be 
generated later should be included into the theory from the very beginning. 
Since the TOE describes the physics down to the classical regime all effective 
coupling constants must appear in the theory. In this manner the coupling 
constant space includes GUT, Standard Model, Nuclear, Condensed Matter, 
Solid State and Atomic Physics parameters, too. Say a QCD quark-gluon 
vertex appears as a multi-particle composite vertex on the level of the TOE. 
We must include all composite vertices in the action what might be needed 
to characterize the physics at lower energies. 

Suppose that we increase the energy of the measurement from few GeV. 
The renormalized trajectory converges to the UV fixed point of QCD only 
in the model computations. In the real world the electro-weak interactions 
which are represented by the non-renromalizable current-current vertices well 
below the Standard Model scale start to deflect the trajectory form the QCD 
UV fixed point. As we further increase the energy the local effective vertices 
explode at the threshold of the weak vector bosons. From now on we see 
the scaling laws of the Standard Model. The renormalized trajectory stays 
in the vicinity of the UV fixed point of the Standard Model (the triviality 
problem of the Higgs sector is now ignored for simplicity) so long as the 
effective vertices generated by the exchanges of a superheavy vector boson 
are weak. By repeting this argument at each intermediate fixed point we find 
that the renormalized trajectory of the TOE visits several fixed points, those 
of the GUT, Standard Model, QCD, QED, Condensed Matter, Solid State 
and Atomic Physics among others, as we move towards the IR. Finally there 
is an IR fixed point. There are further fixed points in the coupling constant 
space which can be approached by the renormalized trajectory when the 
environmental parameters, such as the temperature or chemical potentials 
are properly tuned. 

Since the same operator algebra is classified in the vicinity of each fixed 
point it may happen that a given operator turns out to be relevant at some 
fixed point and irrelevant at others. It is certainly right that the renormalized 
trajectory is influenced only by the relevant operators of a fixed point in 
the scaling regime of the fixed point in question. Thus we have islands of 
universality around each fixed point. But a coupling constant which happens 
to be irrelevant at one fixed point may turn out to important at another fixed 



26 



point. Thus it is a quite involved question that what operators are relevant 
at a certain energy range because the answer depends not only on the energy 
range considered but all scaling regime between the energy range in question 
and the true UV fixed point. We shall see below the possibility that other 
fixed points toward the IR direction may influence the result, too. 

3.3 From the superconductor to the Standard Model 

To simplify the mixing of different fixed points consider the strong and 
the electromagnetic interactions for electrons, muons and nucleons at finite 
baryon density. This gives a well defined UV scaling regime until the Landau 
pole of QED. For the appropriate choice of the chemical potential we may 
have another scaling regime at lower energy which is related to Condensed 
Matter or Solid State Physics. The two scaling regimes define two classifica- 
tion schemes for the operators. Let us call an operator relevant or irrelevant 
at a fixed point if it is contained in the relevant scaling operator set of the 
fixed point or not, respectively. Each operator belongs to one of the classes 
(rel,rel), (rel,irr), (irr,rel) or (irr,irr). Here the first and the second property 
refers to the behavior of the operator in the UV or the IR scaling regime. 
The electron mass, m e , is of the type (rel,rel) since it is a renormalizable pa- 
rameter of the QED lagrangian and influences the phase transitions at lower 
energies. On the contrary, the muon mass, m M , is (rel,irr) because the effects 
of the muons are shielded by those of the electrons at energies below m^. 

The interesting class is the (irr,rel). A coupling constant of that type 
drops as we lower the observational energy from the UV cut-off and becomes 
undetectable. But it starts to grow as we arrive at the IR scaling regime 
and may play important role there and at lower energies. This is a "hidden 
parameter" in the sense that it must be specified at the microscopical scales 
but its effects appear at much longer length scales only. It represents an 
elementary interaction which is not detectable on the level of the elementary 
constituents and becomes important only when certain long range structure 
such as the solid state lattice is formed. An example of this type of operator is 
an effective four fermion contact term which incorporates the weak attractive 
forces between the electrons of the solid due to the phonon exchange. This 
non-renormalizable term is responsible for the superconducting ground state 
whose influence becomes dominant at very low energy. This four fermion 
contact term generated by the phonons is the analogy of the contribution 
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0(J 2 ) in (H). 

According to the strategy of the Wilson-Kadanoff blocking all operators 
which are generated at lower energy scales should be included in the hamilto- 
nian from the very beginning. The appearance of non-renormalizable terms 
should not be a serious problem except for the UV fixed point of the TOE. In 
fact, since all lower energy theory is effective only the scaling laws inferred in 
the vicinity of a fixed point change as we go up in energy and the apparent 
explosion of the non-renormalizable coupling constants gives rise to stability 
in the inter-fixed point region. One is left with wild speculations only con- 
cerning the TOE. A possible scenario to include all coupling constant is to 
require the UV finiteness of this ultimate theory. 

The crucial question whose clarification requires the detailed quantitative 
analysis is whether the initial value for the coupling constants of the type 
(irr,rel) effect the infrared behavior of the theory. In other words, whether 
the dynamical growth of these coupling constants as we pass the crossover 
in lowering the cut-off between the two fixed points is modified by the initial 
conditions taken at the UV side. In case of QED this is the question whether 
the Fermi constant or other effective vertex of similar structure generated 
at the level of the Standard Model or beyond influences the supercurrent 
density in solids. 

We find an aspect of the condensate formation in the example above which 
is different from those observed in regard to the mini-instantons. The four 
fermion interaction becomes important at energies which are well below the 
characteristic mass scale, m e . What kind of degrees of freedom correspond 
to this new operator ? One part of the answer is trivial. Namely, the degrees 
of freedom which are responsible of this vertex belong to the coherent state 
of photons with nonzero momentum, the solid state lattice. The excitations 
of this condensate, the phonons, generate the effective vertex. The other 
part of the question, the degrees of freedom influenced by this vertex is less 
obvious. They come from a condensate again, the boose condensate of the 
Cooper pairs. Thus we have two condensates below m e and their dynamics 
are related to a non-renormalizable operator in the lagrangian. Other field 
theoretical models show similar behavior |29| . 



Historically, the discovery of the superconducting state of matter was a 
total surprise. The understanding of Quantum Field Theory was not suffi- 
cient at that time to predict such phenomena by a theory which was well 
tested in few particle process only. The renormalization group provides us 
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the language and the framework to study problems like this and the sys- 
tematic search of relevant but non-renormalizable operators can reveal such 
surprises. 



3.3.1 Wegner-Haughton equation 

In order to separate the impact of different scaling regimes on the renormal- 
ized trajectory one has to follow the mixing of each operator which might 
prove to be relevant at some fixed point or important at the crossowers. The 
differential form of the renormalization group equation, PS |, is more suited 



for this goal because it handles the mixing of infinitely many coupling con- 
stants in a very economical manner. We shall derive the equation for a scalar 
field theory by using sharp momentum space cut-off. 

Denote the bare action by Sk[<f>] where k is the UV cut-off. We shall 
obtain Sk~vk[(p] in the loop expansion, 

e S k -vM = f £)[0] e -g»frNj = e -^[^o]-|trIn ^Mol +0(Vfc/fc) ; ^ 



where the Fourier amplitude of the fields (f){x) and <p(x) is non- vanishing for 
p < k — Vfc and k — Vk < p < k, respectively and the saddle point is given 
by 5S ^^ = 0. Since the n-loop contributions include an n-fold integration 

over the functional space they are proportional to (6k/k) n in the absence of 
massless singularities in the given kinematical region. Thus 5k /k appears as a 
new small parameter and the exact functional differential equation obtained 
in the limit VA; — > includes the one-loop contribution only. 
In order to simplify fl65|) we use the gradient expansion, 



oo „ 

<?[</>] = E dxUni^x)^ 2 * 1 ), (66) 

n=0 J 

where U n is a homogeneous function of order 2n in the derivative. This 
expansion will be truncated at n — 1, 

S[4>] = Jdx^-id^ + U^), (67) 

and furthermore the simplification Z(<f)) = 1 will be used to derive a simple 
differential equation for the potential U. In order to pick up the local poten- 
tial from the action we choose a homogeneous background field <ft(x) = $. 
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The functional differential equation reduces to 



e -Vl/ fc _ v *(*) = e -^*(*)-¥/^ln(P a +n"(*)) > 



(68) 



where the F stands for the space-time volume and the integration extends 
over the shell k — Vk < p < k. In the limit VA; — > one easily finds 



where Qa denotes the d-dimensional solid angle. This equation represents 
the one-loop resummed mixing of the coupling constants of the potential 
Uk($) = Y,n9n& n - In fact, the expansion of the logarithm in the second 
derivative of the potential gives 



up to a field independent constant. This is the sum over the Feynman graphs 
contributions which come from the infinitesimal loop integration volume. The 
circumstance that the right hand side includes the running potential 
rather than the bare one, U\(<&), indicates that the contributions of the 
successive eliminations of the degrees of freedom are piled up during the in- 
tegration of the differential equation and the solution of the renormalization 
group equation resummes the perturbation series. The solution of the differ- 
ential equation interpolates between the bare and the effective potential as 
k is lowered from the original cut-off A to zero. 

3.4 IR fixed point 

The infrared fixed point is always trivial for theories with mass gap. In fact, 
as the block size extends beyond the correlation length the evolution of the 
coupling constants slows down and we find a manifold of stable fixed points. 
The only relevant coupling constants is the mass since it is divergent in the 
units of the cut-off. But this controls a quadratic operator so has trivial 
effects only. In order to find an example where the IR fixed point generates 
non-renormalizable relevant operators we look for models with massless ex- 
citations. The massless one component 4 model is not appropriate due to 
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(69) 




(70) 
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the Coleman- Weinberg mechanism We need a symmetry to keep the 

mass gap zero in the presence of the interactions. The simplest example is 
the linear sigma model in the symmetry broken phase. Let us start with the 
O(N) invariant lagrangian, 

L = \{d^ a f + Um\ (71) 

a = 1, • • • , N, in four dimensions. The renormalization group equation for 
the potential is of the form |]31] 

kd k U = - (ln(A; 2 + d 2 e U) + (N - 1) ln(A; 2 + d 2 tr U)^ , (72) 

where dt and d tr denote the derivatives in the along the vacuum expectation 
value and the transverse directions of the internal space, respectively. The 
beta function of the n-order massive mode vertex is defined by 

8 n = kd k d%U k (<S>). (73) 

Note that these functions depend explicitly on the scale k and $. The 
strength of the effective interactions for the particlelike excitations are ob- 
tained by setting $ a =< <f) a (x) >. 

Let us identify the leading IR piece of the beta functions. Since 

8 fe 4 ( w*m | (N x d&u*m \ (7A) 

Pl " i67r 2 U 2 + a|^(<i>) + l 'k 2 + diu k {^))' { ' 

the most important IR contribution of the higher order beta functions comes 
from the highest order power of the transverse denominator, 

«■ = ( - 1 »" (Ar - ^(toIm)"* 1 + oVie&D). (75) 

The dimension of the corresponding coupling constant, g n , is 4 — n so the 
beta function, /3 n , for g n = g n k n ~ A in the IR regime is 

«• = '-^wC tSw/t )"' 1 + °^ S ^U)) + (» - Dfc. (76) 

According to the Goldstone theorem 9% r U k= o(< <p >) = 0. We shall assume 
that dndj r U k (^) > and dj r U k (< 4> >) = o(k) which is supported by the 
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simple one-loop solution. The result is that the odd vertices whose beta 
function is negative correspond to relevant coupling constants. The vertices 
<p n , n = 6, 8, • • • which are irrelevant at the UV scaling regime give rise 
relevant operator(s) at the IR fixed point of the symmetry broken theory. 
The inclusion of the renormalization of Z(<$>) leaves our conclusion unchanged 

& 

3.5 Couplings of the Goldstone modes 

The IR divergence in the beta functions poses an interesting problem for the 
chiral perturbation expansion and the Standard Model. It is well known that 
the effective theory of the Goldstone modes, the nonlinear sigma model is 
IR finite. This comes about because the Goldstone modes fluctuations have 
no restoring force by symmetry so they interact with each other via gradient 
couplings which suppress the IR divergences. This makes the chiral models 
which include only the Goldstone modes IR finite. 

It seems surprising at the first moment that the IR finiteness is lost when 
the massive modes are added to the Goldstone particles. Though the on- 
shell amplitudes remain IR finite off-shell the IR divergences appear in the 
massive particle Green functions. These divergences can be found in the one- 
loop effective potential of the linear sigma model where the heavy particle 
legs are connected by a massless particle loop. This is just the contribution 
singled out in fl75|). This divergence comes from the coupling between the 
Goldstone and the heavy modes via the potential U((p a ). This coupling con- 
tains no derivatives and leads to the IR divergence in the beta function. The 
absence of the derivatives in this coupling can be understood by noting that 
the transverse modes of the linear sigma model are not exactly the Goldstone 
modes but contain the heavy mode at higher order in the fluctuations around 
the vacuum. Since the massless propagator of Goldstone modes comes from 
the second order contributions of the action and the mixing between the 
Goldstone and the heavy modes is of higher order the massive mode inter- 
acts by itself by the long range massless propagator. This IR divergence of 
the effective potential makes the dynamics of the condensate sensitive for 
the coupling constants which would have been completely unimportant oth- 
erwise. Similar strongly coupled long range interactions have been noticed 
in the framework of our extended Nambu-Jona-Lasinio model, after equation 
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In regard to the Standard Model this argument raises the possibility that 
the higher order, non-renormalizable Higgs vertices which are generated by 
the exchange of the superheavy particles of the GUT scale may have an un- 
usual strong influence on the low energy physics. These vertices are certainly 
small as we follow the renormalized trajectory around the energy scale of the 
Higgs mass, M#, but they may become large in the IR scaling regime. Since 
we have a differential equation the initial conditions at the UV cut-off deter- 
mine the effective coupling constants in the IR region. Both the suppression 
in the UV and the amplification in the IR generate qualitatively similar 
power dependence around the fixed points according to the linearized block- 
ing equations. So it seems plausible to expect seizable sensitivity of the IR 
physics on these non-renromalizable bare parameters. If the detailed numer- 
ical solution of the renormalization group equation reveals such a sensitivity 
then new parameters of the Standard Model are found. These parameters 
become important well below the mass gap. What degrees of freedom do 
they influence ? The answer points again to the condensate since < <f>(x) > 
is determined by the long range interactions between the heavy and the Gold- 
stone modes. It is reasonable to assume that the number of new parameters 
is the number of the dynamically generated condensate. There is a differ- 
ence with superconductivity in the manner the non-renormalizable vertices 
influence the excitation spectrum. The four fermion interaction leads to the 
Cooper pair formation whose condensate gives rise a new low energy excita- 
tion spectrum. In the case of the Higgs particle < 4>(x) > generates the mass 
for the fermions and the gauge bosons. Thus the modification of the strength 
of the condensate feeds back to the whole mass spectrum of the theory. 

4 Conclusions 

The confinement of quarks has generated different models and led to the 
development of a number of non-perturbative mechanisms. The present lec- 
tures intend to introduce few new pieces to this collection. These are the 
breakdown of the fundamental group symmetry, the role of the nontrivial 
integration measure in the path integral, the presence of singular configura- 
tions in renormalized theories and the appearance of new parameters by the 
condensates or the multiple fixed point structure. QCD has a lesson for us to 
learn in each of these directions. Such a rich and many-sided theory can only 



33 



raise our determination to aim at synthesis and arrive at a comprehensive 
understanding. 
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